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ABSTRACT 

We use oblate axisymmetric dynamical models including dark halos to determine the 
orbital structure of intermediate mass to massive Coma early-type galaxies. We find a 
large variety of orbital compositions. Averaged over all sample galaxies the unordered 
stellar kinetic energy in the azimuthal and the radial direction are of the same order, 
but they can differ by up to 40 percent in individual systems. In contrast, both for 
rotating and non-rotating galaxies the vertical kinetic energy is on average smaller 
than in the other two directions. This implies that even most of the rotating ellipticals 
are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral 
axisymmetric toy models we show that flattening by stellar anisotropy maximises 
the entropy for a given density distribution. Collisionless disk merger remnants are 
radially anisotropic. The apparent lack of strong radial anisotropy in observed early- 
type galaxies implies that they may not have formed from mergers of disks unless the 
influence of dissipational processes was significant. 

Key words: galaxies: kinematics and dynamics - galaxies: elliptical and lenticular, 
cD - galaxies: formation 



1 INTRODUCTION 

The way in which a galaxy has assembled its stars is 
reflected in the distribution of stellar orbits. For ex- 
ample, collisionless Ai'-body collapse simulations predict 
a p redominance of radial orbits in the final remnant 
(|van Albada 1982). In contrast, collisionless galaxy merger 
simulations predict a variety of o rbital compo s itions, de- 
pending on progenitor properties jBarned 19921: Hernguistl 
1992, 199^), the merging geometry llWeil fc Hernauistlll99^ 
Dubinski 199 ^, the progenitor m ass ratios ( Naab fc BurkertI 
2003 ; J esseit. Naab fc BurkertI 2005') and the presence of 
dissipational com ponents (|Naab. Jesseit fc BurkertI I2OO6I : 
I Jesseit et al.ll2007h . Stars in galaxies are approximately col- 
lisionless and the orbital structure - once a galaxy has ap- 
proached a quasi-steady-state - is conserved for a long time. 



To a certain extent then, the assembly mechanism of early- 
type galaxies can be constrained from their present-day or- 
bital structure. 

A global characteristic of the distribution of stellar or- 
bits is its anisotropy. Traditionally, anisotropics of elliptical 
galaxies have been inferred from the {v/cr, e) diagram. In par- 
ticular, the rotation of bright ellipticals has been shown t o 
be insufficient to account for their flattening (|Binnevlll97"8l ). 
However, whether fainter, fast-rotating ellipticals are flat- 
tened by rotation is less easy to determine from the {v/a, e) 
diagram, because isotropic as well as anisotropic systems can 
rotate. In fact, fully general axisymmetric dynamical mod- 
els recently have revealed an anisotropic orbital structure 
in ev en the flattest, fast rotating objects l|Cappellari et all 
l2007f ). One goal of this paper is to investigate numerically 
the connection between anisotropy, rotation and flattening 
in spheroidal stellar systems. 
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In addition, we present global anisotropics for a sam- 
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pie of Coma early-type galaxies. These anisotropies are de- 
rived by analysing long-slit stellar absorption line kinemat- 
ics with axisymmetric orbit models. Our dynamical mod- 
els include dark matter halos. Previous anisotropy deter- 
minations for larger samples of ellipticals (including dark 
matter) were restricted to round a nd non-rotating sys- 
tems, assuming spherical sym metry (|Gerhard et al.lT200ll : 
iMagorrian fc Ballantvnj|2001^ . Spherical models do not ac- 
count for galaxy flattening. In the simplest case, a flattened 
system is axially symmetric. Early axisymmetric models, 
however, did not cover all p ossible anisotro pics (and orbital 
structures, respectively; e.g. ICaroUo et afl 1995). Fully gen- 
eral, orbit-based axisymmetric dynamical models have so far 
only been applied to the inner regions of ellipticals and the 
orbital analysis was made under the a ssu mption that mass 
follows light (e.g. lGebhardt et al.ll2003l and lCappellari et al.l 
[2OO6). By the mass-anisotropy degeneracy, the neglect of 
dark matter could translate to a systematic bias in the corre - 
sponding orbital structure (e.g. lForestell fc Gebhard3l2008h . 
Comparison of anisotropies derived with and without dark 
matter will allow one to quantify such a possible bias. 

We also discuss anisotropies derived from modelling 
mock observations of synthetic A'^-body merger remnants. 
One motivation to do so is that dynamical models of ax- 
isymmetric systems may not be unique. For example, the 
deprojection of an axisymmetric galaxy is intrinsically de- 
generate (jRvbickl ,1987.). Uncertainties in the intrinsic shape 
thereby propagate int o uncertainties on the derived masses 
and anisotropies (e.g. iThomas et al.|[2007ah . Moreover, the 
reconstruction of an axisymmetric orbital system is sus- 
pected to be further degenera te with the recovered mas s 
(e.g. the discussion in Valluri. Merritt fc EmsellemI |2004) . 
The case for a generic degeneracy, beyond the effects of 
noise and incompleteness of the data, is still uncertain (e.g. 
[Magorrian 20o3). Numerical studies of a few idealised ax- 
isymmetric toy models indicate degeneracie s to be moderate 
when modelling realistical ly noisy data sets l|Krainovic et al.l 
l2005l . lTh omas et al.ll2005l ). Since we know the true structure 
of our A^-body modelling targets, we can extend on these 
studies and further investigate potential systematics in the 
models over a broader sample of test cases. 

Another motivation to model A'^-body merger rem- 
nants is to probe whether ellipticals have formed by merg- 
ing. This requires a comparison of the orbital structure in 
real ellipticals with pr e dictions of A^-body s imulations (e.g. 
iBurkert fc NaabI l2005l : iBurkert et all |2008| ). However, be- 
cause of the symmetry assumptions in models of real galax- 
ies, it is not straight forward to compare intrinsic proper- 
ties of A'^-body simulations with models of real galaxies. To 
avoid the related systematics, we here compare models of 
real galaxies with similar models of synthetic A'-body merger 
simulations and both are indicative for true differences be- 
tween real galaxies and merger predictions. 

The galaxy and A''-body merger samples and the mod- 
elling technique are briefly outlined in Sec. [2] Toy models 
of various flattening and anisotropy are discussed in Sec. |3l 
The anisotropies of real galaxies are presented in Sec. U and 
compared with models of A'^-body merger remnants in Sec. [5] 
Implications for the formation process of early-type galax- 
ies are discussed in Sec. [H] and we summarise our results 
in Sec. [7] The influence of regularisation and the inclusion 
of dark matter halos on reconstructed galaxy anisotropies 
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Table 1. Summary of COMA galaxy anisotropics. ( 1-2): g alaxy 
id (GMP numbers from [Godwin, Metcalfe & PcachI Il983r) : (3): 
intrinsic ellipticity e; (4-6): anisotropy parameters 5, f} and 
7 (cf. equations I5I7I I of the best-fit dynamical model; (7): 
(ti/cr)*, i.e. (umax/cTo) normalised by the approximate value 
\/ fobs/ (1 ~ ^obs) of a-u (edge-on) isotropic rotator with the same 
flattening. Note that (v/cr)* is an observable, i.e. i t combines ob- 
served ellipticities t^bs (from column 10 of Tab. 1 of lMehlert et al] 
[2000) and observed velocities i^max and ctq, without reference to 
any dynamical model. 

is discussed in App. [X] In App. [B] we briefly discuss the 
connection between anisotropy and the shape of the circular 
velocity curve in maximum entropy models. We assume that 
the Coma cluster is at a distance of 100 Mpc. 



2 DATA AND BASIC DEFINITIONS 

A complete description of a stellar system is given by its 
distribution function / (DF; the density in 6-dimensional 
phase-space). In a steady-state system the DF / depends 
on the phase-space coordinates only through the integrals of 
motion (jLyriden-BoU 1962). Axisymmetric potentials, which 
are considered here, admit the two classical integrals of mo- 
tion energy (E) and z-component of the angular momentum 
(Lz). In addition, many orbits in astrophysically relevant po- 
tentials are characterised by another, non-classical, so-called 
third integral (I-^ : .Contopouloai 1963f ) . Since integrals of mo- 
tion label orbits and vice-versa, a steady-state system can 
be viewed as a superposition of orbits, each with constant 
phase-space density. Let f{i) denote the phase-space den- 
sity along orbit i, then the total amount of light w{i) on the 
orbit equals 'w{i) — f{i) x V{i) {V{i) is the orbit's phase- 
space volume). The DF - or the weights of a suitable orbit 
superposition model - determine the spatial density p and 
intrinsic velocity dispersions a via 

P = J fd'v (1) 
and 



3 



/ {vi - Vi){vj - Vj) d^v, 



(2) 



(3) 



In the following we will only consider i,j € {z, x, <j), R}, 
where z is the short axis of the density distribution, </> is 
the azimuth around this axis, a; is a fixed Cartesian coordi- 
nate parallel to the equatorial plane and i? is a cylindrical 
radius. Let 



(4) 



denote the totaQ (unordered) kinetic energy in coordinate 
direction i, then the global anisotropy of an axisymmetric 
stellar system can be quantified, for example, by the ratios 
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/3 



and 



(6) 



(7) 



ijCappellari et al1l2007l ). In axisymmetric systems, the three 
anisotropy parameters are related via 

2/3-7 



S = 



7 



(8) 



Non-rotating, isotropic spherical systems as well as classical 
isotropic rotators obey 5 = P — j — 0. 

The DF of real galaxies is not known, but has to be 
reconstructed from photometric and kinematic observations. 
In the next two subsections we will describe the two samples 
of real and simulated galaxies discussed in this paper and 
will briefly outline our modelling method. 

2.1 Coma early- types 

Our sample of observed galaxies (C OMA in the fo l lowing) 
consists of 19 Coma early- types from lXhomas et al.l (|2007bl . 
I2OO8I ). It comprises 2 central cD galaxies, 10 ordinary gi- 
ant ellipticals and 7 SO or intermediate galaxies with lu- 
minosities — 20.3 ^ Mb ^ —22.56 (a single fainter object 
with Mb ~ —18.8 is also included in the sample). The 
galaxies are drawn from the luminosity limited sample of 
iMehlert et al] l|2000l ) and are distributed all over the cluster. 
High-resolution radial profiles of surface brightness, elliptic- 
ity and isophotal shape p arameters 04 an d ae (up to ai2 in 
some cases; cf. iBender fc MoUenhofjl 19871 for a definition of 
the isophotal shape parameters) derived from a combination 
of HST and ground-based imaging were used to calculate 
the deprojected 3d luminosity distribution for several incli- 
nations. The photometric data are complemented by long- 
slit stellar absorption line kinematics along 2 — 4 position- 
angles per galaxy. The kinematic data consists of radial 



profiles of mean velocity, velocity dispersion and higher- 
order moments of the line-of-sight velocity distribution and 
reach out to 1 — 4 reft . Details about the p hotometric and 
kinematic dat a can b e found in j0rgensen fc Franx, (19941). 
iMehlert et"al] 1I2OOOI). IWegner et all (|2002l ). ICorsini et all 
(l2008l ) and lThomas et al.l (|2008D . 

These data were modelled with our implementation 
of Schwarzschild's (1979) o rbit superposition technique 
for axisymmetric po t ential s ([Richstone fc Tremaind Il988l : 
iGebhardt et al.ll2000l . l2003l : iThomas et all 120041 '). Fbr each 
galaxy, we probed for a variety of mass models, composed 
of a stellar mass density (from the deprojected light profile) 
and a parametric dark halo profile. The parameter space 
for the mass models spans the inclination, the stellar mass- 
to-light ratio and the dark halo parameters. In each trial 
potential the best-fit orbit model is calculated by maximis- 



S - ax 



(9) 



where quantifies deviations between observed and mod- 
elled kinematicfl The function 



S = -2_^w{i) In^)^ 



w(i) 

n(i) 



(10) 



is used to smooth the orbit models. In the absence of 
any other constraints the maxim isation of S y ields orbital 
weights w{i) cx Q{i) (.Richsto ne fc Tremaind ll988), such 
that the yet not specified can be regarded as weight- 
factors for the w{i). When modelling real galaxies or mock 
observations of Ai'-body merger remnants, we assume that 
there is no preferred region is phase-space and each orbit 
is given an a priori-weight equal to its phase-space volume: 
= V{i). Then, 



5" = -^m(i)ln 



- / /ln/d\d^u 



(11) 



equals the Boltzmann entropy, which drives models towards 
a constant density in phase-space. 

The (binned) deprojected luminosity density is used as 
a boundary condition to solve equation (|9} and the regu- 
larisation parameter a in equation ( O has been calibrate d 
by means of Monte-Carlo simulations (jThomas et al.|[2005l ). 
The final, best-fit orbit model is obtained from a x'^-analysis. 



2.2 Af-body merger remnants 

We have applied the same modelling code to mock observa- 
tions of synthetic A-body merger simulations. In brief, we 
have modelled six merger remnants, each projected along its 
three principal axes (models of projections along the long, 
intermediate and short axis of the merger remnants will be 
shortly referred to as X, Y and Z-models later on). The six 
merger remnants are taken f rom the sample of colli sionless 
disk-l-bulge-l-halo mergers of iNaab fc BurkertI ()2003l ). They 
have mass ratios between 1:1 and 4:1 and sample the en- 
tire distribution of intrinsic shapes and orbital make-ups, 
including extreme cases. An orbital analysis of the A-body 



^ In the following we will only consider anisotropics (and ki- 
netic energies, respectively) inside the effective radius, because 
this is the radius inside which kinematical observations are typi- 
cally available to constrain the orbital structure of real galaxies. 



^ We use the full information of (binned) line-of-sight velocity 
distributions when fitting real galaxies and A^-body merger rem- 
nants. 
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systems is given in ljesseit. Naab fc Burker3 (|2005l ). We have 
simulated typical Coma observations for each projection: the 
merger remnants were placed at a distance of 100 Mpc and 
photometric and kinematic profiles with similar resolution 
and spatial coverage as in the COMA sample have been ex- 
tracted. For a detailed discussion o f the models the reader 
is referred to iThomas et al.l l|2007al ). 



3 TOY MODELS 

Oblate stellar systems can owe their shapes to a variety 
of difi'erent orbital configurations. Classically, one has of- 
ten distinguished between two proto-typical cases: flatten- 
ing by rotation and flattening by anisotropy. Thereby, flat- 
tening by rotation is used to term an otherwise round and 
isotropic system which appears flattened - and rotating - by 
extra-light on near-equatorial, high angular-momentum or- 
bits (populated with the same sense of rotation). Flattening 
by anisotropy refers to systems with a depression of stars 
with high velocities perpendicular to the equatorial plane 
(Jizz < Hflfl ~ n^^). However, in fact there are infinitely 
many orbit superpositions that accoun t for a given galaxy 
shape. Some of these are discussed in iDehnen fc GerhardI 
l| 19931 ). Different orbital structures can be distinguished by 
their different anisotropics. In the following we will numeri- 
cally construct (self-consistent) toy models that are designed 
to (1) reproduce a given, flattened, density distribution ex- 
actly, but (2) have different intrinsic anisotropics. 



3.1 Self-consistent models without rotation 

The models are orbit-based and similar to those described 
in Sec. [2] However, here we only require the models to re- 
produce a given density (a = in equation [9]) . Various ex- 
pressions for the factors Q{i) in equation pop will be used 
to impose different anisotropy structures (see below). 

For our simple toy models we assume a stellar density 



p oc m ^(m -I- 1) ^ 
jHernquisdlfggol ) with 



2 r.2 , Z 



(12) 



(13) 



Equations (|12|) and (|13p describe systems with constant flat- 
tening q. They approximate the light profiles of elliptical 
galaxies reasonably well. 

Flattening and maximum entropy. Let fs denote the 
DF that maximises the entropy of equation ((TTJ subject 
to the density constraints. The squares in the top panel 
of Fig. [1] illustrate the connection between anisotropy and 
flattening for fs'- the three panels show the anisotropy pa- 
rameters from equations (|5I7|I as a function of the intrinsic 
ellipticity e = 1 — q (cf. equation [13}. While S and /3 increase 
with flattening, 7 is roughly constant. In maximum-entropy 
models the flattening thus arises from a suppression of en- 
ergy in z-direction, while the balance between the energies 
in R and <p is roughly conserved. In this sense, the maximum 
entropy models fs resemble the classical case of flattening 
by anisotropy. The only difference is that 7 7^ (cf. App. [B] 



for a discussion of 7) . Note that we calculated the toy mod- 
els with the same library setup as used for the Coma galaxy 
models. 

Flattening by a classical f — f{E, L^). A classical two- 
integral DF f2i, which only depends on E and Lz, can be 
approximated via equations (|9I10|I with 



n{i) = 



1 - C{i) 



C{i) = V{i) I ^ V{j) 



(14) 



(15) 



and 

J{i) = {je{l,...,N}:Lz{j)^Lz{i), EU)^E{i)} (16) 



(A'' is the total number of orbits) . Equation (I14p derives from 
the constraint that for f2i ~ f{E, Lz), the phase-space den- 
sity of any orbit i with energy E{i) and angular momentum 
Lz (i) has to equal the mean phase-space density of all orbits 
j with the same E{i) = E{j) and Lz{i) — Lz{j), i.e. 



w{i) X] 



V{i) 



(17) 



This case is included in Fig. [T]by the circles. That the 
r2(j) from equation lfT4|l indeed yield / « fiE, Lz) is demon- 
strated by /3 « 0. The flattening of the corresponding sys- 
tems comes from an excess energy in 0-direction with respect 
to the isotropic case (7 < 0; orbits with high angular mo- 
mentum are strongly populated). The relationship between 
5 and e is similar as in maximum entropy models. 

Note that DFs / « f{E,Lz) develop noticeable phase- 
space density peaks on orb its with high angular momentum 
(|Dehnen fc GerhardI Il994h . It is likely this property that 
lowers their entropy as compared to the fs models. Flat- 
tening by anisotropy mainly involves shell orbits which ap- 
proach closely the intrinsic minor-axis. Their phase-space 
volumes are much larger than those of equatorial near- 
circular orbits with high angular momentum. Even a small 
change in the phase-space density along shell orbits can re- 
duce the amount of light near the minor-axis considerably 
and, thus, result in a significant flattening. The larger frac- 
tion of phase-space involved in this type of flattening, com- 
pared with a strong overpopulation of the relatively small 
region in phase-space occupied by near-circular orbits (as in 
cases where / = f{E,Lz)) explains why objects which are 
flattened by anisotropy have the higher entropy. 

Flattening with radial anisotropy. Model DFs /^ad ob- 
tained with 



^{i) = ['"apo(i) - rpcri{i)f X V{i) 



(18) 



are biased towards orbits with a large difference rapo — f peri 
between apocentre and pericentre radius (radially extended 
orbits). Such models are radially anisotropic (/3,7 > 0; cf. 
triangles in Fig. [l| . The relationship between 5 and e is again 
similar as in the previous models. 

The latter is no surprise, as for self-consistent ellipsoids 
with constant flatt ening, 6(e) can be calculated from the ten - 
sor virial theorem (|Robertdl 19621 : iBinnev fc Tremaindll987l ) : 
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Figure 1. Anisotropy parameters 5, /3 and 7 (as labelled in the panels) versus intrinsic ellipticity e for the toy models /g (squares), 
/27 (circles) and /^ad (triangles). Toy models of the same type are connected by dotted lines. Top: models without rotation; bottom: 
models with rotation (/.t = 0.9; cf. equation I24I I: light symbols: without dark matter halo; solid lines: tensor virial theorem applied 
to oblate spheroids with constant flattening. 



where 



and 



0.5 



i(e) — e\J\ — 



e(l — e^)~^-^ — arcsin(e) 



(19) 



(20) 



e = ^l-il-eY. (21) 

The solid line in the upper-left panel of Fig.[T]shows relation 
l|19p . Our numerically constructed orbit models follow this 
line well. 

Note that, if A'' DFs fi project each to the same spatial 
density, then any convex linear combination f\ = Xj fj 
with '^Xi — 1 will do so. The properties of f\ will be inter- 
mediate between those of the individual fi . 

3.2 Rotation 

The just discussed toy models (and any linear combination 
of them) are non-rotating, because in our choices for f2(i) we 
haven't distinguished between prograde and retrograde or- 
bits. A large variety of rotation patterns can be constructed 
from any DF / as follows: each orbit in an axisymmetric 
potential comes in two flavours, one prograde (with posi- 
tive Lz > 0) and one retrograde (Lz < 0). Both share the 



same spatial shape but differ only in the sign of the velocity 
component around the axis of symmetry. Thus, the spatial 
density will only depend on the sum 



f+{E,Lz,h) = f{E,Lz,h)+f{E,-Lz,h 



(22) 



of hght on corresponding prograde and retrograde orbits. 
The amount of rotation, instead, will depend on the differ- 
ence between the population of the prograde and retrograde 
orbits, respectively. This can be quantified, for example, by 
the fraction /j, of light on the prograde of each orbit pair: 



fiE,Lz,h) 



(23) 



(Lz ^ 0). For simplicity, let's assume from now on that /i is 
the same for all orbits. Then, any 



f^{E,Lz,l3 



fif+{E,Lz,l3) ■■ >0 

(l-^)/+(iJ,-Lz,/3) : Lz<0 



(24) 



with n G [0, 1] (to remain positive definite) will give rise to 
the same density profile as / (/^ = /+), but with different 
degrees of internal rotation. For example, in case of = 
1/2 prograde and retrograde orbits are populated equally 
and there will be no rotation in the corresponding system. 
With fi = 1 {fi = 0) only prograde (retrograde) orbits are 
populated (maximum rotation). 
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Figure 2. First three panels, from left to right: elliptieity tobs versus classical vja for the toy models /2/ (grey circles), /s (grey 
squares) and /jad (gr^y triangles). For each toy model the cases = 0.5, 0.7, 1.0 (no rotation, intermediate and maximum rotation) 
are shown. Models with the same /i are connected by the dotted lines. All toy models are edge-on such that = e. Outer right 
panel: e^ha versus v/a for Coma galaxies (filled circles: edge-on galaxies, open circles: non edge-on galaxies). For comparison, the 
maximally rotating f2i and /g models are also shown in the outer right panel. In all panels, the dashed lines approximate edge-on 
isotropic rotators (cf. equation 125( 1. 



The bottom panel of Fig. \T\ shows anisotropies for the 
toy models of Sec. I3.1l with fi = 0.9. While (3 is independent 
of the amount of rotation, S decreases and 7 increases with 
increasing rotation. The latter reflects that in our toy models 
the total energy in ^-direction is constant. Any increase of 
the ordered motion is thus to the expense of a smaller £70. 

Fig. [2] illustrates where the toy models appear in the 
{v/a, e) diagram. The figure shows the three cases p — 0.5 
(no rotation), /i = 0.7 (intermediate rotation) and fj, = 1.0 
(maximum rotation). On the y-axis, the ratio Wmax/cro of 
the maximum rotation velocity (vmax, along the projected 
major-axis) and the central velocity dispersion (ao, aver- 
aged inside rcs/2) is shown. All models are edge-on. The 
highest rotation rates at a given flattening are obtained 
for f2i, because of its strongly populated high angular mo- 
mentum orbits (n^^ > TIrr « H^^). However, fs models, 
which are not flattened by an excess of light on high an- 
gular momentum orbits (relative to the isotropic case) but 
instead by a suppression of orbits with large ^-velocities 
(IIzz < Urr « Tl^g>) can reach (wmax/cro) « 0.75 as weU. 
The dashed lines in Fig. [2] approximate classical isotropic 
rotators by 

(^) (25) 

l|Kormendvlll982l ). Up to e ~ 0.4, fs models can appear 
in the same region as classical isotropic rotators, although 
they are not flattened by rotation in the classical sense (e.g. 
/3 7^ 0) . Radially anisotropic models are dominated by orbits 
with low angular momentum and have generally low rotation 
rates. 

A complete picture of an axisymmetric galaxy's flatten- 
ing mechanism requires knowledge of the amount of rotation 
(e.g. v/a) and at least one anisotropy p arameter (e.g . /3, 7 or 
6 or the parameter a in the notation of lBinnevll2005l ). Alter- 
natively, two anisotropy parameters also specify the global 
orbital structure. In any case, the full information about the 
anisotropy and the flattening mechanism cannot be provided 



by the {v/a, e) diagram alone. For example, four among the 
five intrinsically most fiattened COMA early-types are very 
close to the isotropic rotator line in Fig. [2] However, they 
are shaped by a combination of /3 > and 7 < (cf. Tab.[T|. 



3.3 Influence of a dark matter halo 

The presence of dark matter around a galaxy affects the 
shape of the stellar orbits. Some of the models just discussed 
may not exist, if an additional dark matter halo reshapes the 
potential significantly. To check this, we have recalculated 
all our toy models in a potential, where a spherical, logarith- 
mic dark halo has been added to the stellar potential. The 
parameters of the halo (its core radius and its asymptotic 
circular velocity) have been set according to the dark mat- 
ter s caling relations in COMA early-types (jThomas et al.l 
|2008| ). The derived anisotropies in the new potential differ 
in no case by more than 0.1 from the original ones (but 
see the discussion in App. [Bj. Especially, the relationship 
between /3 and e, that arises from the maximisation of the 
orbital entropy also appears in potentials with a realistic 
dark matter halo. This does not necessarily imply that the 
neglect of dark matter in models of real galaxies has no ef- 
fect on the derived anisotropies, because it may enforce a 
redistribution of the orbits (cf. next Sec. [Jjl . 



4 REAL GALAXIES 

Fig. [3] shows the connection between anisotropy and flat- 
tening in real galaxies. The intrinsic flattening of COMA 
galaxies is expressed in terms of 

f dR SB(R) e(R) , , 

JdRR^SBiR) 

(|Binnevl[2005l ). Here, R is the radius along the projected 
major-axis and SB(i?) and e{R) are the surface-brightness 
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Figure 3. Top: Anisotropy parameters S, (3 and 7 (as labelled in the panels) versus intrinsic ellipticity e for real galaxies (cf. mid 
panel in top row). Bottom: same for models of A'^-body merger remnants. Light: toy models from equation l|27|l . connected by dotted 
lines. Solid lines in left panels: tensor virial theorem applied to oblate spheroids with constant flattening. 



profile and ellipticity profile in the edge-on projection. For 
an axisymmetric system (with flattening q) e « (1 — 5). 
Lines in Fig. [3] trace three different toy models 

/a = A/5 + (l-A)/2, (27) 

(cf. SecO the three models are designed to rotate by using 
^ = 0.75 in equation [24)1 . DFs / « f[E, L^) are inconsistent 
with the global orbital structure of most galaxies (because 
(9 > in observed galaxies) . Most galaxies have orbital prop- 
erties between those of fs and f2i (with some rotation). 

Fig. [3] also includes a nisotropi es and flattenings of 24 
early- types from Cappclla ri et al.l (|2007). These galaxies 
are a subsample of the 48 Es/SOs of the SAURON sur- 
vey (|de Zeeuw et al.ll2002l ) , which uniformly covers the plane 
of observe d flattening Coba and M b (for Mb < —18). The 
galaxies of lCappellari et al] (|2007l ) are drawn from this sur- 
vey according to various requirements, among them con- 
sistency with axial symm etry (according to 2d kinemati- 
cal maps). The galaxies of ICappellari et al] l|2007l ') (shortly 
SAURON in the following) are on average fainter than the 
COMA galaxies. 

Although the samples do not match exactly, the 
anisotropies of COMA and SAURON galaxies are found in 
the same range. However, the COMA sample contains rel- 
atively more anisotropic but nearly round galaxies on the 
one hand and more highly flattened but isotropic galaxies 



{5 ~ 0) on the other. As a result, the trend for 5 and /3 
to increase with e which is seen in the SAURON sample is 
not obvious when considering the complete COMA sample 
(even not if the two Coma galaxies with the most uncertain 
anisotropies are ignored - the two central galaxies GMP2921 
and GMP3329). 

The relation between {3 and e is weaker in the COMA 
galaxies in part due to a few round but anisotropic galax- 
ies - for example GMP1750 and GMP5568 with e ^ 0.2 
and 13 ~ 0.26 — 0.36. Both ga laxies show weak minor-axis 
rotation (|Thomas et al. I l2007bh and could be slightly triax- 
ial systems. In addition to differences among nearly round 
galaxies, anisotropies in COMA and SAURON galaxies also 
slightly differ at high e. The latter is most clearly seen in 
5 versus e: two highly flattened COMA galaxies (GMP1990 
and GMP2440, e « 0.6) have 5 ^ Q. One of these galaxies is 
likely close to edge-on, because of its high ob served elliptic- 
ity (m ax Sobs ~ 0.625, cf. the radial proflle in lThomas et al.l 
l2007bl ) and its signiflcant isophotal shape distortions. We ex- 
pect the model of GMP1990 to be well constrained, because 
of the near edge-on inclination (minimal uncertainties in 
the deprojection) and its far-extendin g multi-slit kinematic 
data. For the other galaxy (GMP2440') iMehlert etHI (|2000l ) 
quote only a modest observed ellipticity tobs = 0.33 at r^s 
and the intrinsic flattening comes mostly from the low incli- 
nation of the model. Note that this galaxy is far above the 
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isotropic rotator line in the right panel of Fig.[2](GMP2440 is 
the only non edge-on galaxy above the isotropic rotator line). 
A maximum-entropy like DF is ruled out for this galaxy, be- 
cause even the maximally rotating version of the fs model 
would not allow for the high observed rotation rate. Thus, 
even if we would have underestimated the inclination of 
this system. Fig. [2] shows that its orbital structure must 
be significantly deviant from maximum-entropy models. All 
in all then, modelling uncertainties are unlikely to explain 
the outstanding anisotropies of G MP1990 and GMP244q . 
In fact, a comparison with Fig. 3 in ICappellari et al.l (|2007| ) 
reveals that the SAURON sample does not include galax- 
ies like GMP1990 and GMP2440, because (1) for only one 
object the observed ellipticity is significantly larger than 
Cobs > 0.5 (NGC4550) and (2) even the fastest rotators in 
the SAURON sample are closer to the isotropic rotator line 
than GMP2440. 

In addition to differences in the sam ple selection also 
the m odelling methods differ in the details. ICappellari et al.] 
l|2007l ) use similar orbit-based dynamical models as we do 
here, but SAURON anisotropies are calculated inside a fixed 
aperture with a radius of 25" . A fixed aperture encloses dif- 
ferent fractions of the stellar mass in different galaxies, de- 
pending on system size and distance. For the COMA galax- 
ies we give anisotropies inside r^s- In some galaxies local 
anisotropies vary significantly with radius ([Thomas et al.l 
l2007bl ). such that the radius of comparison is crucial. In ad- 
dition, SAURON models are based on the assumption that 
mass follows light. As it has been stated already in the in- 
troduction, the assumption of a constant m ass-to-light ra- 
tio can result in artificially large ^-energies (|Thomas et al.l 
l2005l . lForestell fc Gebhardtll2008l ) or low 7, respectively. Re- 
garding Fig. [31 SAURON galaxies do not have systemat- 
ically lower 7 than COMA object s. For the only two ex - 
ceptions (NGC4473 and NGC4550) ICappellari et all (|2007l ) 
report evidence for counter-rotating, disk-like components 
that likely cause their large ^-energies. The small effect that 
the neglect of dark matter has on the anisotropies likely 
reflects the fact that we only consider anisotropies aver- 
aged inside r < VeS, where the assu mption that mass fol- 
lows light is most clo sely fulfilled (e.g. iGerhard et al.ll200ll . 

omas et al.ll2007bl ). For the Coma galaxies a quantitative 
comparison of models with and without dark matter is made 
in App.lXl 

The spatial coverage with kinematic data in the in- 
ner regions is sparse in the COMA galaxies (long-slit data) 
compared to the SAURON objects (2d kinematical maps). 
In regions of phase-space that are not well constrained by 
the observed kinematics, the dynamical models are mainly 
driven by regularisation. Thus, because the spatial coverage 
is lower in the COMA galaxies, their anisotropies could be 
biased towards the adopted regularisation scheme. Specifi- 
cally, COMA g alaxy models are reg ularised towards max- 
imum entropy l|Thomas et al.|[2005^ . However, the middle 
panel in the top row of Fig. [3] does not show any bias of the 
COMA models towards the maximum entropy relation. In 
fact, SAURON galaxies are on average closer to this relation 
than COMA galaxies. This indicates that regularisation is 
not the main driver for the COMA galaxy models. Also, in 
App. we give an explicit comparison of COMA galaxy 
models with standard and with weak regularisation. We do 
not find significant differences. 



Both the intrinsic ellipticity and the anisotropy depend 
on the inclination of the models. For the COMA galax- 
ies, we probe three different inclina tions and use the one 
that fits best ((Thomas et al.|[2007lj ). while inclinations for 
SA URON galaxies are deriv ed from two-integral Jeans mod- 
els (|Cappellari et all |2006| ). The inclination is best con- 
strained for highly flattened galaxies, because these have 
to be close to edge-on. For three of the COMA galaxies 
(GMP0756, GMP1176 and GMP1990) large ellipticities to- 
gether with significantly discy/boxy isophotes indeed indi- 
cate close to edge-o n inclinations (for example, GMP1176 
exhibits > 10; ICorsini et al.l |2008| ). In contrast, two 
among the five galaxies with e ^ 0.5 owe their flattening 
in part from the relatively low inclination of the best-fit 
model (GMP0282, GMP2440; cf. Tab. (U. These galaxies 
provide the smallest and largest anisotropies, respectively, 
at high e (cf. middle panel in the top row of Fig. [3]). This 
suggests that the method to determine the inclination for 
the COMA galaxies does not result in a specific bias of the 
derived anisotropies. 

We conclude that slight differences between the 
SAURON and the COMA anisotropies are mostly due to 
the different sample selections, while differences in the mod- 
elling methodology (including differences in the data cover- 
age) seem to be negligible. 



5 COMPARISON OF OBSERVED GALAXIES 
AND iV-BODY MERGER REMNANTS 

The lower panels of Fig. [3] display the models of TV-body 
merger remnants (cf. Sec. O). In terms of 5 versus e and (5 
versus e, these models do not differ strong ly from models of 
real galaxies (see also lBurkert et ahlbOOSl ). However, while 
7 in models of merger remnants, 7 is often negative in 
models of real galaxies. Is this discrepancy in 7 indicative 
for the merger remnants having a different orbital struc- 
ture than real galaxies, or does it merely reflect systematics 
caused by the symmetry assumptions in our models? 

Reconstructed and true intrinsic anisotropic^ and flat- 
tenings of the merger remnants are compared in Fig. 
The one merger remnant closest to oblate axial symmetry 
(OBLATE), is reconstructed with high accuracy from the X 
and Y-projections (edge-on). This is plausible, because for 
this remnant the assumption of axial symmetry is a good 
approximation. Furthermore, in the edge-on case the depro- 
jection becomes unique and the intrinsic degeneracies in the 
dynamics are likely smallest. 

However, the general trend in the axisymmetric models 
is to underestimate both, the flattening and the anisotropy 
of the merger remnants. X and Y-projections allow a better 
reconstruction of shape and anisotrop y than Z -projcctions. 
It has already been discussed in Th omas et al.l (|2007a l that 
the assumption of axial symmetry enforces an inclination 
mismatch in the Z-models: while the triaxial remnants ap- 
pear flattened in the Z projection (face-on), axisymmetric 



Note that while Ilxx = ^yy in axisymmetric systems, Hxx ^ 
in the merger remnants. For the intrinsic 5 of the merger 
remnants we use the average ijlxx + nyy)/2 instead of Hxx in 
equation ((Sj. 
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Figure 4. From top to bottom: anisotropy parameters <5, /3 and 7 versus (intrinsic) ellipticity e. Heavy symbols: intrinsic parameters 
of Af-body merger re mnants (abbreviation s refer to the TRIAX, PROLATE, ROUND, FLAT, ELONG and OBLATE merger 
remnants discussed in iThomas et~alll2007ah : light: Schwarzschild models of merger remnant projections. From left to right: models 
of X, Y and Z-projections are highlighted in colour. In each panel, dashed lines indicate which model belongs to which merger 
remnant. Solid lines in top row: (5(e) of equation II19I I. as in Fig. [T] 



systems are necessarily round when seen face-on. Then, be- 
cause the models are forced towards a wrong viewing- angle, 
(1) the intrinsic flattening is underestimated and (2) X, Y 
and Z-axes of models and remnants do no longer correspond 
to each other. For example, a Z-mode l's 7 measures a differ - 
ent energy ratio as the remnant's 7 (|Thomas et al.|[2007al ). 
Had we compared the Z-models with the apparent shape 
of the remnant in Z-projection and with the energy ratios 
along axes of models and remnants that correspond to each 
other, then the differences would have been much smaller 
(for example Ae < AS < 0.1). 

Real galaxies are seen at random viewing angles. Start- 



ing from our models of principal projections it is difficult 
to predict directly the analogous distributions of 5{e) or 
/3(e) for the realistic case of random projections. However, 
because the projections along principal axes yield extreme 
kinem atical and photometrical propert ies of the merger rem- 
nants (|jesseit. Naab fc BurkertI |2005| ). it can be expected 
that dynamical models of projections along intermediate 
viewing-angles will have properties intermediate between 
those of the models from principal projections. We have ver- 
ified this for two out of the six merger remnants (ELONG 
and OBLATE) by modelling additional 11 projections (at 
intermediate viewing-angles). Assuming that this result can 
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Figure 5. Correlations among tlie anisotropy parameters 5, (3 and 7. Top: models of real galaxies (squares, triangles and circles: toy 
models from equation 1271 connected by dotted lines); bottom: models of merger remnants (squares, triangles and circles separate 
models of projections along different principal axes as indicated in the lower-left panel; crosses with circles: intrinsic anisotropics 
of the Af-body remnants). 



be generalised to other remnants as well, then Figs. [3] 
and 2] suggest the following: if real galaxies would resem- 
ble the modelled merger remnants, then one would see ap- 
proximately the same relationships 5(e) and /9(e) as in the 
COMA and SAURON galaxies. However, it is clear from 
Fig. |4] that 7^0, for a sample of randomly projected ob- 
jects like our modelled A^'-body merger remnants. Thus, in 
this respect, models of many real galaxies differ from our 
comparison sample of synthetic A'^-body merger remnants: 
models of merger remnants are always radially anisotropic 
(n_Rij > ~ Tlzz), while models of real galaxies are char- 
acterised by n0<^ > Urr > Uzz. 

This fact is further illustrated in Fig. (5] which shows 
correlations among the anisotropy parameters. Even though 
shape and anisotropy cannot be recovered simultaneously 
(in some cases), the anisotropy correlations in the models of 
the merger remnants and in the merger remnants itself are 
very similar to each other. Again, the main difference be- 
tween real galaxies and merger models is the offset between 
their 7 distributions. 

Besides the fact that merger remnants have on average 



positive 7 > 0, while real galaxies have 7 « (on aver- 
age). Fig. [S] shows that the distribution of anisotropics in 
the merger remnants is tighter than in real galaxies. This 
may reflect the similarity in the initial conditions of the 
Ai'-body simulations (most noteworthy the similarity in the 
progenitors and the fact that we only consider coUisionless 
mergers) . 



6 IMPLICATIONS FOR THE FORMATION 
PROCESS OF EARLY-TYPE GALAXIES 

The anisotropy parameters defined in Sec. |3] are only global 
measures of the orbital structure. A full understanding of 
the formation process of early-type galaxies can only be pro- 
vided by spatially resolved anisotropy profiles. For example, 
equatorial near circular orbits obey, in the epicycle approx- 
imation, the local relation 

_ 1 / dhrUcirc \ 
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wher e Wcirc is the circular velocity l|Binnev fc Tremaind 
1 19871 ). In a typical galaxy potential the circular velocity 
curve is flat (dlnucirc/dlnr ~ 0) and equation (|28|) pre- 



dicts cr^ ~ 0.5 crj. Since the epicycle approximation holds 
for perturbed rotating disks, we do not expect the majority 
of early-types in our sample to be well described by equation 
l|28p . However, it might be relevant for the most flattened, 
rotating and discy objects in our sample. Instead, at least 
some of these (for example GMP1176 and GMP3958) have 
negative 7 (i.e. > an). This does not rule out a disk heat- 
ing scenario for these galaxies, however, because locally we 
find al ^ (0.5- •• 0.7) near the equatorial plane in these 
ga laxies (around r m r eg; cf. the radial anisotropy profiles 
in lThomas et al.ll2007bl ). 

In case of the coUisionless A'^-body merger simulations, 
already the averaged anisotropy parameters reveal signifi- 
cant differences to the models of real galaxies. Which phys- 
ical processes are responsible for this discrepancy? 

The orbital structure of the models of merger remnants 
is largely driven by a population of cen tral box orbits in the 
Af-body systems ijThomas et al.llioOTal 'l . They cause the cen- 
tres of the merger remnants to become triaxial/prolate and 
are, for example, largely responsible for the wrong viewing 
angle of the Z-models. Dissipation during a merger can have 
a significant effect on t he shape and the projected properties 



of the final remnant ( Barnes & Hcrnquist 1996; Cox et a,l.l 
20061: iRobertson et ahfe oOG: Naab. Jcsscit & Burkcrt 2003; 
Jesseit et al.l 12007 ). Already 10 percent of gas are suffi- 



cient to suppress central box orbits and to produce an 
approximately axisymmetric remnant in binary mergers 
l|Naab. Jesseit fc Burkerj|2006l ') (but this result is based on 
simulations without star formation). 

Multiple, simultaneous minor mergers likewise produce 
remna nts less triaxial than co Uisionless binary merger rem- 
nants (|Weil fc Hernauis3ll996l ), but the corresponding kine- 
matics have not yet been studied in detail. Successive minor 
merging does not necessarily lead to different final remnants, 
at least if the cumulative mer ged mass becomes similar 
to th e most massive progenitor i Bournaud. Jog fc Combea 
I2OO7I ). Again, detailed predictions for the orbital make-up 
and the shapes of the line-of-sight velocity distributions have 
not yet been worked out. 

Note that the central dark matter densiti es in COMA el- 
liptic als are larger than in present-day spirals (jThomas et al.l 
|2008| ). Even if ellipticals have formed by some variant 
of merging, present-day spiral galaxies are unlikely the 
progenitors for the bulk of giant el liptic als (see also 
iNaab fc Ostriked 120071 ). iBurkert et all (|2008l ) pointed out 
that A'^-body systems, which have assembled h ierarchically 
in their cosmological simulations (|Naab et al.|[2007, ), or by 
binary mergers with star-formation and black-hole feedback 
are consistent with the trend between S and e in observed 
galaxies. 



7 SUMMARY 

We have discussed the relationship between anisotropy and 
flattening in toy models, in models of real galaxies, in merger 
remnants and in models of merger remnants. Models of ob- 
served galaxies generally exhibit /3 > and 7 « 0. We do 



not find strong correlations of the anisotropy parameters 5, 
P and 7 with intrinsic ellipticity e. 

In toy models with maximum entropy for a given den- 
sity distribution we find /3 to increase with e, while 7 < 0. 
Observed galaxies appear close to these maximum-entropy 
relations, but exhibit a large degree of individuality. Rota- 
tion appears in anisotropic (/? > 0) as well as isotropic sys- 
tems (/3 ~ 0), suggesting that the flattening of the galaxies 
largely arises from a suppression of stars with large ener- 
gies perpendicular to the equatorial plane. This is similar to 
the classical notion of flattening by anisotropy and rules out 
DFs / ~ f{E,Lz) for most early-type galaxies. 

The global similarity between models of observed galax- 
ies and our maximum-entropy toy models suggests that 
early-type galaxies are largely relaxed stellar systems. How- 
ever, there are differences in the details that probably con- 
tain valuable information about the assembly mechanism of 
the galaxies and will be addressed in a future paper. 

Numerical simulations indicate that both strongly ra- 
dially anisotropic (7 — » 1) and strongly tangentially 
aniso tropic systems (7 ^ 0) can become unstable 
(e.g. iMerritt fc StiaveUil Il990l : ISellwood fc MerrittI Il994l : 
iNipoti et al.1 |2002| ). Maximum entropy models have in- 
termediate anisotropics and are likely stable. Thus, the 
anisotropics of observed galaxies may not only be under- 
stood as being the most likely ones (in the sense of yielding 
the maximum entropy at a given flattening) but could also 
reflect stability constraints. So far we lack detailed studies 
exploring the stability of axisymmetric systems with dark 
matter halos and various intrinsic anisotropics. Since our 
(three integral) toy models can be easily transformed to 
A'-body systems (cf. I Thomas et al.l l2007ah they provide a 
suitable tool to setup both artificially anisotropic as well as 
realistic and observationally motivated initial conditions. 

In models of real galaxies the unordered kinetic en- 
ergy in the azimuthal direction, H^^, can exceed the ra- 
dial energy Hrr by up to 40 percent. This separates real 
galaxy models from similar models of coUisionless A'^-body 
binary disk mergers, which are instead characterised by ra- 
dial anisotropy {Hrr > H^^ ^ lizz)- Because we have ap- 
plied the same modelling machinery to both, the real galax- 
ies as well as the synthetic A'-body merger remnants, our 
findings indicate a true difference between their intrinsic 
properties. Especially, we have shown that if real galaxies 
would resemble our merger remnants, then corresponding 
dynamical models of real data would be radially anisotropic, 
irrespective of the systematics introduced by the assumption 
of axial symmetry. 

The radial anisotropy of the merger remnants is related 
to a population of central box orbits. Because dissipation 
during a merger can efficiently suppress box orbits, our re- 
sults suggest that dissipation played an important role dur- 
ing the formation of intermediate mass to massive early-type 
galaxies. 

In this paper we focussed on the comparison of real 
galaxies with coUisionless binary disk merger simulations. A 
similar analysis, but for gaseous mergers with star formation 
and/or for galaxies formed in cosmological simulations could 
give more insight into the actual formation paths of elliptical 
galaxies. 
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APPENDIX A: THE INFLUENCE OF 
REGULARISATION AND DARK MATTER ON 
RECONSTRUCTED GALAXY ANISOTROPIES 

The Coma galaxy models discussed in this paper are de- 
rived using our standard regularisation parameter a = 0.02 
(cf. equation |9]). The strength of the standard regularisa- 
tion has been calibrated by means of Monte-Ca rlo simula- 
tions of isotropic rotators ([Thomas et al.l[2005h . To check 
how much the choice of a affects our results, we (1) deter- 
mined the best-fit dynamical model at q = 1 and (2) re- 
calculated the anisotropies of all galaxies from these weakly 
regularised models (at a = 1, the minimum is usually 
reached). In the top row of Fig. IAU we show both models 
with standard and with weak regularisation for comparison. 
As can be seen, lowering the regularisation has almost no 
effect on the derived anisotropies. Especially, there are still 
at least four galaxies with significantly negative 7 < 0. 

In the bottom row of Fig. lAll we make a similar com- 
parison for models with and without dark matter: squares 
indicate the anisotropies of our best-fit models with a con- 
stant mass-to-light ratio (no dark matter halo). As expected, 
when assuming that mass- follows- light, 7 become smaller 
(the amount of iji-energy is increased to compensate for the 
missing dark mass). From the bottom- right panel one would 
expect that the average 7 becomes negative when the ra- 
dial increase of the mass-to-light ratio (caused by a dark 
halo) is neglected. This is not the case in the SAURON 
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Figure Al. As Fig. [Sj but in the top row Coma galaxy models are plotted with standard regularisation (filled circles; a = 0.02 
in equation lOj and with weak regularisation (stars; a = 1). The bottom row shows Coma galaxy models with (filled circles) and 
without (squares) dark matter halos (in both cases the standard regularisation was applied). 



sample, however, although ICappellari et aP l|2007h assumed 
that the mass-to-light ratio is constant with radius in their 
models. That neglecting dark matter has a stronger effect 
in Coma galaxies is likely related to the fact that our kine- 
matical data reach out into the region where dark matter 
becomes noticeable (> rgff), which is probably not the case 
in many SAURON galaxies (where the data extend only out 
to < roff). 



APPENDIX B: THE RADIALLY RESOLVED 
ANISOTROPY STRUCTURE OF 
MAXIMUM-ENTROPY TOY MODELS 

The maximum-entropy toy models fs discussed in Sec. 13.11 
resemble the classical flattening by anisotropy, except that 
they are only approximately isotropic in R and To investi- 
gate where this anisotropy comes from, we have constructed 
maximum-entropy toy models in potentials that include a 
dark matter halo. The halo density distribution is assumed 
to follow 



jDehnenlllQQ^ ). To mimick reahstic halos (cf.lTh omas et al.l 
l2007bl ) we choose a flat central density core (r] = 0.05) and 
we set the flattening q of the halo equal to the flattening of 
the luminous component of the toy model (cf. equation [13]). 
We investigated three mass models: (1) no halo, (2) a mass 
model that has a roughly constant circular velocity curve 
and (3) a mass model with a rising Vdic in the outer parts 
of the model. The corresponding circular velocity curves for 
q = 0.7 are shown in the bottom panels of Fig. IBll 

The upper panels of Fig. IBll display the radial profiles 
of the local anisotropies 



/3cq(r) 



and 



7cq(r-) = 1 - 



o-fl(r)2 



MO! 



(B2) 



(B3) 



Pdm oc m ''(m -I- 1)'' 



(Bl) 



along the equatorial plane (averaged within ^ ±11.5°, 
where -& is the latitude). As one can see, the anisotropy in 
the meridional plane (/3cq) does not depend on the shape 
of the gravitational potential. Thus, the relation between (3 
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Figure Bl. Radial anisotropy profiles /3cq (top) and 7cq (mid- 
dle) for maximum-entropy models (/g; flattening q = 0.7) with 
dark matter halos. Left-hand side/squares: halo tuned to re- 
sult in an approximately flat circular velocity curve; right-hand 
side/triangles: halo leading to an increasing circular velocity in 
the outer parts of the model. For comparison, the case without 
halo is shown on both sides (circles). Circular velocity curves of 
the models (scaled to the maximum circular velocity vq without 
halo) are shown in the bottom panels. 



and e is largely independent from the gravitational potential 
and closely related to the entropy maximisation. 

However, beyond rcir, where dark matter starts to influ- 
ence the shape of the circular velocity curve, 7oq is different 
in the three different potentials. Fig. IB2l shows that the local 
value of 7eq - along the equatorial plane - is directly con- 
nected to the logarithmic slope Ovcirc of the circular velocity 
curve. In general then, because 7 from equation ((Tjl is the 
spatial average of 7oq (and the local anisotropics along other 
position angles in the meridional plane), its exact value is 
not set uniquely by the entropy maximisation but also de- 
pends on the shape of the circular velocity curve. In practice, 
however, deviations with respect to the model without halo 
become noticeable only beyond reft, such that even the spa- 
tially averaged 7 of the toy models does not depend strongly 
on whether a halo is included or not. 

Note that the relation revealed by Fig. IB2I is different 
from the epicycle relation (|28p . This is expected, because the 
azimuthal velocity dispersion in the toy models largely 
results from the fact that they do not rotate. Instead, the 
dispersion predicted by the epicycle approximation arises 
from perturbations on circular orbits in a rotating disk. 
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Figure B2. Local anisotropy 7oq (along the equatorial plane) 
versus logarithmic slope a^irc = din tidrc/d Inr of the circular 
velocity curve. Symbols and colours as in Fig. IBll 



